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ABSTRACT 


Radiation pattern prediction analysis of elliptically polarized cavity 
backed aperture antennas in a finite ground plane is carried out using a combined Finite 
Element Method (FEM) /Method of Moments (MoM)/ Geometrical Theory of Diffraction 
(GTD) technique. The magnetic current on the cavity-backed aperture in an infinite 
ground plane is calculated using the combined FEM /MoM analysis. GTD, including the 
slope diffraction contribution, is used to calculate the diffracted fields due to both soft 
and hard polarizations at the edges of the finite ground plane. Numerical results for the 
radiation patterns of a cavity backed circular spiral microstrip patch antenna excited by a 
coaxial probe in a finite rectangular ground plane are computed and compared with 
experimental results. 


2 



LIST OF SYMBOLS 


v 


Vr 




P cA’ P cB 


Mi 

CO 

D 

S 

D, 


h 


E 


H, 


H 


H 

i 

K 

H 1 


ap 

inp 

e// 

Wf 

d 

0A 

d 

(f)A 


Kb 


H' 


C| )fi 


H , 


ec 


Del operator 
Relative permittivity 
Relative permeability 
Permeabilty of free space 
Intrinsic impedance of free space 
Defined in equation (15) 

As defined in equations (16) and (17) 

Angular frequency 
Soft diffraction coefficient 
Hard diffraction coefficient 
Electric field 

Electric far field, 0 component 

Electric far field, (f> component 

Magnetic field at the aperture 

Magnetic field at the input plane 

Magnetic far field, 0 component 

Magnetic far field, (f> component 

Diffracted magnetic field at point A, 0 component 

Diffracted magnetic field at point A, (f> component 

Diffracted magnetic field at point B, 0 component 

Diffracted magnetic field at point B, (f> component 

Diffracted magnetic field at point C, 0 component 


3 



H, 

H‘ 


dfC 
d 

0 D 

d 

<S)D 

H ‘vf 

i 


M 


M x 

M 

y 

M x 



h 


r, 0, <]) 

s ap 

e 

inp 


y 

T 


V 


W v w 2 

X,Y,Z 


Diffracted magnetic field at point C, (|) component 

Diffracted magnetic field at point D, 0 component 

Diffracted magnetic field at point D, (f> component 

Total magnetic far field after adding the diffracted field, 0 component 

Total magnetic far field after adding the diffracted field, (f> component 

7=1 


Free space wavenumber 
Fourier transform variable with respect to x 
Fourier transform variable with respect to y 
Magnetic current in the aperture 
x-component of magnetic current 
y-component of magnetic current 
Fourier transform of M x (equation (7)) 
Fourier transform of M y (equation (8)) 
unit normal 
Spherical coordinates 
Return loss(dB) 

Integration over the aperture surface 
Integration over the input plane 
Testing function 
Integration over the volume 
Dimensions of the finite ground plane 
Cartesian coordinates 


4 



1. INTRODUCTION 


Cavity backed antennas are important in aerospace and automobile 
applications due to their conformal nature. Previously a combined FEM/MoM/GTD 
technique has been applied to linearly polarized cavity backed antennas in a finite 
ground plane[l]. In this method, the far field due to the magnetic current in the aperture 
is considered to be incident on the edge of the finite ground plane to produce the 
diffracted field. In GTD, when the electric field is normal to the edge of the ground plane, 
it is referred to as hard polarization and when it is tangential to the edge of the ground 
plane it is referred to as soft polarization. For soft polarization, the electric field incident at 
the diffraction point on the ground plane is zero. In such cases, GTD may be used with 
slope diffraction. Slope diffraction is based not on the magnitude of the incident field at 
the point of diffraction, but rather on the slope magnitude of the incident field at the 
point of diffraction. 

The application of slope diffraction to cavity backed antennas is given in 
detail in reference^]. In the present report, a combined FEM/MoM method is used to 
analyze the elliptically polarized cavity backed antennas in a finite ground plane. GTD, 
with slope diffraction coefficients is used to augment the radiation patterns of these 
antennas on a finite ground plane. A particular case of a cavity-backed Archimedean 
spiral microstrip antenna fed by the center conductor of a coaxial line is considered. For 
this elliptically polarized antenna, GTD along with slope diffraction is used to account 
for the radiation patterns of both E Q and E () polarizations. The radiation patterns of the 


5 



spiral antenna are calculated and compared with experimental results measured at 
NASA Langley Research Center's Low Frequency Antenna Measurement Facility Good 
agreement between the calculated and measured data is obtained. 
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2. THEORY 


The geometry of the problem to be analyzed is shown in figure 1. For a 
linear, isotropic and source free region the electric field satisfies the vector wave equation 


Vx| -VxlY 
Mr ) 


2 

W 


(i) 


where (I and £ are the relative permeability and permittivity respectively of the 
medium in the cavity The time variation exp (jcot) is assumed and suppressed 
throughout this paper. 


To facilitate the solution of the partial differential equation in (1) for the 
cavity via FEM, equation (1) is multiplied by a vector testing function T and the result 
is integrated over the volume of the cavity By applying suitable vector identities, 
equation (1) can be written in its weak form[8] as 
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Applying the divergence theorem to the integral on the right hand side of 
equation (2) and using Maxwell's equation VxE = -j(Si\x o \xH , where u , is the 
permeability of free space and co is the angular frequency, equation (2) can be rewritten 
as 

III (3) 

V Mr y s 


7 



where h is the unit outward normal to the surface and H is the magnetic field on the 
surface of the cavity. As the tangential electric field is zero on the perfect electric 
conducting (PEC) walls of the cavity, the surface integral in equation(3) is nonzero only 
over the aperture opening in the infinite ground plane and the input aperture. 

Ill -( v x^J*(vxlJdv-^e r ||jf«£dv-ya)^J \[Txh)*H ap ds 

V» r V s ap 

= j'cop o J J [rx hj • Hjnpds (4) 

s, np 

where 5 is the aperture surface and S u is the input surface (see fig. 1). H ap i s the 
magnetic field at the aperture and Hi np is the magnetic field at the input surface. 

At this point, the problem can be divided into three parts, the first part 
involving the discretization and evaluation of volume integrals on the left hand side of 
equation (4) and the second part involving the discretization and evaluation of the 
surface integral over S . The third part involves the surface integral over the input 
aperture. The volume integral is evaluated by discretizing the volume of the cavity with 
tetrahedral elements and evaluating the integral over the volume of each tetrahedral 
element. These element volume integral contributions are added together over all 
tetrahedral elements to form a sparse matrix. The second integral is evaluated over the 
surface of the aperture. The discretization of the volume of the cavity into tetrahedral 
elements automatically discretizes the aperture surface into triangular elements. 
Assuming the current basis functions similar to that described in [3], an integral 



equation is formed, which contributes to the dense part of the system matrix. Image 
theory is used to account for the presence of the infinite ground plane. The third integral 
is evaluated over the input coaxial aperture surface and forms the excitation column 
matrix. Solving for the system matrix formed by these three integrals, the electric fields 
in the cavity are obtained and hence the magnetic current M in the aperture S . 

Using the magnetic current M ( x , y ) over the aperture, the magnetic far 
field is given by the following equations: 

jk exp (-jkr) 

H qjj ( r, 0, (>) = — cosQ\_M x (k x ,k y ) costy + M y (k x ,k y ) sin§] (5) 

jk exp {-jk r) _ 

0. ♦> = -- 2 „ \M, (*.-• V CO s* -M, (k x , k y ) sin*] (6) 

where 

M x (k x , k y ) = J J M x (x, y) exp {j {k x x + k y y) } dxdy (7) 

and 

M y [k x , k y ) = J J M y (x, v) exp {j {k x x + k y y) } dxdy (8) 

k = ^ O cos(|)sin0 (9) 

k = sin 4> sin 0 (10) 

Using equations (5) and (6), radiation patterns in the forward direction of a 
cavity backed antenna in infinite ground plane are computed. The effect of the finite 
ground plane is incorporated by employing the edge diffraction coefficients along with 
the slope diffraction coefficients. The radiation fields are supplemented by the additional 


9 



field contributions arising from the diffraction by the edges of the finite ground plane. 
The diffracted field supplements the field in the forward direction, whereas it accounts 
totally for the field in the shadow region (behind the finite ground plane) where the 
infinite ground plane solution gives a zero field. 

The geometry of the finite ground plane is illustrated in figure 1. The 
midpoints on the two edges of length W 2 are given by A and B, respectively, whereas the 
midpoints on the two edges of length Wj are given by C and D, respectively. The total 
diffracted field at A and B are given by [2,6] 


<4 ~ l °) D S (PcA> V, P', n) +J i 


dH, 


e // 


30 


d ’lV 


M_s 

3\|/' 


pcA.^.y.P' 0 .«y 


P cA 
r 


1/2 


exp (- jk r ) exp (Jk ( d sin0) 


( 11 ) 


Kb ~ \ l n)D s (p cB , \ 2t V|A P' 0 . n) + i- 


})H 


e // 


30 


. 7t 

d ’2’*J 


Ms 

3\|/' 


PcB- %2> V. P' D - n) 


P cA 
r 


1/2 


exp (- jk r ) exp (- jk sin 0) 


(12) 


for (f> polarization, and 
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for 0 polarization. In the above equations. 
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71 

P' o = - for normal incidence[6, pp.811] 


n=wedge factor[6, figure 13-30] 


Similarly, the diffracted fields at C and D are also calculated by replacing 
p cC = p cD = d = — . The detailed expressions for diffraction coefficients D g/h are 
given by [6, eq.13-89]. The slope diffraction coefficients are given in [2]. The total 
diffracted field is given by 
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( 18 ) 


and 



Ka + h \ 
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(|) = 0° plane 
(f> = 90° plane 


j H Ia + Kb ^ = 0° plane 
^ H'' ¥: + H^ n § = 90° plane 


(19) 


The diffracted fields for (f> = 180° and (f> = 210° planes can also be 
calculated in a similar manner. The total radiated far field is given by adding the 
diffracted magnetic field and magnetic far field given in equations (5) and (6); i.e., 


H 'vf= H »ff+K ( 2 °) 

(21) 

The radiated electric far field is given by 

£ e = nflV (22) 

E, » = 11 H' w (23) 


Using the magnetic current in the aperture in equations (5) and (6) and evaluating 
equations (20),(21), (22) and (23), the radiation patterns of the cavity-backed aperture 
antennas with a finite rectangular ground plane are computed. 
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3. NUMERICAL AND EXPERIMENTAL RESULTS 


To verify the validity of the analysis presented, a cavity-backed spiral 
microstrip antenna in a finite rectangular ground plane is considered. The spiral 
microstrip antenna is designed to produce elliptical polarization. The geometry of the 
spiral antenna is shown in figure 2. A two and a half turn, two arm circular 
Archimedean spiral is designed to be fed at the center by the center conductor of the 
coaxial probe. The substrate is 1/16" thick with a relative dielectric constant of 2.33. The 
width of the spiral is designed to be 2/23". The microstrip spiral antenna is placed in a 
2" diameter cavity of 1/16" depth with a rectangular ground plane of 24" length and 12" 
width. 

3.1 Experimental Procedure: 

The spiral antenna was fed at its center with the center conductor of 
RG-141 coaxial semi-rigid cable. The copper shield was conductively bonded to the 
copper clad back side of the antenna. This back side provided a conductive path to the 
ground plane when placed in the machined cavity in the ground plane. The antenna/ 
ground plane assembly was then placed on the roll-over-azimuth-over-elevation 
positioning system in the NASA Langley Research Center Low Frequency Antenna 
Measurement Facility. 

Electrical boresight was accomplished by sighting the characteristic null, 
normal to the plane of the antenna, of a spiral fed in the manner described. Calibration of 
the system was not performed as the data were normalized. 
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Principal plane measurements were accomplished for both polarizations, 
horizontal E Q ) and vertical) T 0 ), and for two orthogonal positions of the antenna, 
(f> = 0°, 180^ and (f> = 90°, 270 r , shown in figure 2. 

3.2 Numerical Modelling: 

The cavity volume with the spiral antenna is modelled using a commercial 
package COSMOS /M[8] and meshed with 8,740 tetrahedral elements which resulted in 
1820 triangular elements in the aperture surface. The radiation patterns for both 
horizontal (E e ) polarization and vertical ( £ ( ,, ) polarization are computed using 
equations (20) and (21) respectively. Figures 3 (a) and (b) are plots of normalized 
radiation patterns for horizontal polarization at (f> = 0°, 1 80^ and (f> = 90°, 210° 
respectively at 4 GHz. It can be seen that the measured data and computed results agree 
very well. Figures 4 (a) and (b) are plots of normalized radiation patterns for vertical 
polarization at (f> = 0°, 1 80° and (f> = 90°, 210° respectively at 4 GHz. Good agreement 
between the computed and experimental results is found in the forward direction of the 
antenna. There is considerable disagreement in the backward direction. Possible causes 
for this disagreement are discussed in the following subsection. Similarly, computed 
results for normalized radiation patterns for both horizontal and vertical polarizations 
are plotted along with the experimental data in figures 5 and 6, respectively, at 5 GHz. 
The disagreement between the computed and measured results in figure 6 are discussed 
in the following subsections. 


14 



3.3 Experimental limitations : 

Two primary contributions should be mentioned as possible sources of 
measurement error. The first is the poor return loss, 5 U , of the spiral antenna used in the 
measurement. The antenna, fed as previously described with no balun or matching 
network, was matched poorly to the 50 Q, measurement system. A Wheeler can test [9] was 
performed and indicated that most of the return loss that was measurable could be 
attributed to other than radiation resistance. Therefore, only small signal levels were 
coupled to the coaxial line feeding the antenna and caused a degradation of overall 
system dynamic range of the measurement set-up. The second inaccuracy was due to the 
placement of the remote mixer and relatively low signal levels mentioned above. 
Standard practice is to place the mixer as close as possible to the antenna under test. The 
mixer for these measurements was placed behind the ground plane, so that losses at the 
desired test frequency could be minimized. As a consequence, the mixer and the cable 
were directly illuminated for angles of 0 > 90 r ’ . Although the cable and mixer are both 
shielded, neither are perfect. It is felt that this coupled with the already small desired 
signal level on the cable, due to mismatch and the position relative to the main lobe of 
the antenna, can explain some of the differences between the measured and computed 
results for angles (9 >90°). Further experimentation possibilities include designing a 
feed network to properly match the antenna and using absorber around the remote 
mixer, if still necessary. 

3.4 Computational limitations: 

Due to the relatively small thickness of the substrate of the spiral antenna, 
the volume of the cavity has to be meshed with a good quality tetrahedral mesh. To 
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account for the rapid variation of the field in the cavity, it is essential that a few layers of 
elements be present to account for the field variation accurately This requires a large 
number of tetrahedral elements and proportionately large number of triangular 
elements in the aperture surface, which results in a large number of non-zero elements in 
the system matrix. As the frequency increases the electrical depth of the cavity increases, 
since the wavelength decreases, which in turn requires smaller tetrahedral elements to 
meet the field sampling requirement. Hence, the computations are limited to 4 and 5GHz 
for the cases investigated here. With proper meshing and availability of computational 
facilities such as parallel computing, this antenna configuration can be accurately 
analyzed at higher frequencies using the present analysis. 

4. CONCLUSIONS 

An analysis of cavity-backed elliptically-polarized aperture antennas in a 
finite rectangular ground plane is carried out using a combined FEM/MoM/GTD 
technique. Slope diffraction coefficients are successfully used to account for the 
diffracted field of soft polarization, whose incident field at the diffraction point is zero. 
The radiation patterns for a circular Archimedean spiral microstrip antenna are 
computed with a finite rectangular ground plane and compared with the experimental 
data. Within the computational and experimental limitations both computed and 
measured radiation patterns show very good agreement validating the analysis 
presented in this report. 
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INPUT 


Figure 1 Geometry of a cavity backed aperture in finite ground plane 
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Ground plane 
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Substrate 


50F1 coaxial feed 


Figure 2 Cavity backed spiral microstrip antenna in a finite ground plane 
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Figure 4 Radiation Patterns of spiral antenna (figure 2) at 4GHz - Vertically ) Polariza- 
tion 

(a) (])=0 o and cf)=180° 

(b) 4>=90° and d=270° 
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